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ABSTRACT 

Gravitational lensing effects arise from the light ray deflection by all of the 
mass distribution along the line of sight. It is then expected that weak lensing 
cluster surveys can provide us true mass-selected cluster samples. With numer- 
ical simulations, we analyze the correspondence between peaks in the lensing 
convergence K-map and dark matter halos. Particularly we emphasize the differ- 
ence between the peak k value expected from a dark matter halo modeled as an 
isolated and spherical one, which exhibits a one-to-one correspondence with the 
halo mass at a given redshift, and that of the associated K-peak from simulations. 
For halos with the same expected k, their corresponding peak signals in the k- 
map present a wide dispersion. At an angular smoothing scale of = 1 arcmin, 
our study shows that for relatively large clusters, the complex mass distribution 
of individual clusters is the main reason for the dispersion. The projection effect 
of uncorrelated structures does not play significant roles. The triaxiality of dark 
matter halos accounts for a large part of the dispersion, especially for the tail 
at high n side. Thus lensing-selected clusters are not really mass-selected. To 
better predict ^-selected cluster abundance for a cosmological model, one has 
to take into account the triaxial mass distribution of dark matter halos. On 
the other hand, for a significant number of clusters, their mass distribution is 
even more complex than that described by the triaxial model. Our analyses find 
that large substructures affect the identification of lensing clusters considerably. 
They could show up as separate peaks in the K-map, and cause a mis-association 
of the whole cluster with a peak resulted only from a large substructure. The 
lower-end dispersion of k is attributed mostly to this substructure effect. For 
6g = 2 arcmin, the projection effect can be significant and contributes to the 
dispersion at both high and low k, ends. 



Subject headings: cosmology: theory — dark matter — galaxy: cluster — general 
— gravitational lensing — large-scale structure of universe 
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1. Introduction 

Because they are directly associated with the mass distribution of the universe, gravi- 
tational lensing effects are powerful probes of spatial structures of the dark matter. Strong 
lensing phenomena, such as multiple images of background quasars and giant arcs, have 
been used to constrain inner mass profiles of lensing galaxies and clusters of galaxies (e.g., 
Gavazzi et al. 2003; Bartelmann & Meneghetti 2004; Ma 2003; Zhang 2004). Weak lensing 
effects, on the other hand, enable us to study mass distributions of clusters of galaxies out 
to large radii (e.g., Bartehiiann & Schneider 2001). Cosmic shears, coherent shape distor- 
tions of background galaxies induced by large-scale structures in the universe, provide us a 
promising means to map out the dark matter distribution of the universe (e.g., Tereno et al. 
2005; Van Waerbeke 2005). 

Of many important studies on lensing effects, the aspects of weak lensing cluster surveys 
attract more and more attention (e.g., Reblinsky & Bartelmann 1999; White et al. 2002; 
Padmanabhan et al. 2003; Hamana et al. 2004; Haiman et al. 2004). Clusters of galaxies 
are the largest virialized structures in the present universe. Their formation and evolution 
are sensitive to cosmologies, and therefore can be used to constrain different cosmological 
parameters, such as as, fim and the equation of state of dark energy, where as is the rms of 
the extrapolated linear density fluctuation smoothed over 8/i~^Mpc, and Jl^ is the present 
matter density in units of the critical density of the universe (e.g., Bahcall & Bode 2003; 
Fan & Chiueh 2001; Fan & Wu 2003; Haiman et al. 2001) . There are different ways 
finding clusters. The optical identification based on the concentration of galaxies suffers 
severe projection effects. X-ray and Sunyaev-Zel'dovich (SZ) effect are associated with the 
intracluster gas, and have been used extensively in cluster studies (e.g., Rosati et al. 2002; 
Carlstrom et al. 2002). However, most of the theoretical studies concern the abundance 
of clusters in terms of their masses (e.g.. Press & Schechter 1974; Sheth & Tormen 1999; 
Jenkins et al. 2001), therefore it is crucial to get rehable relations between different survey 
observables and clusters' mass. The properties of intracluster gas are affected significantly 
by gas physics, which we have not fully understood yet. Thus there are large uncertainties 
in relating X-ray and SZ effect with the total mass of a cluster. On the other hand, lensing 
effects of a cluster are determined fully by its mass distribution, and therefore clean mass- 
selected cluster samples are expected from weak lensing cluster surveys. 

However, weak lensing surveys have their own complications. Lensing effects are asso- 
ciated with the mass distribution between somxes and observers, and thus the lensing signal 
of a cluster can be contaminated by other structures along the line of sight. The intrinsic 
ellipticities of source galaxies can pollute the lensing map and lower the efficiency of cluster 
detections considerably. Besides and more intrinsically, clusters themselves generally have 
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complex mass distributions, and their lensing effects can be affected by different factors in 
addition to the total mass. 

Therefore for extracting cosmological information from a sample of lensing-selected clus- 
ters, three main theoretical issues need to be carefully studied. Firstly the lensing effects 
from clusters must be understood thoroughly. Secondly the significance of projection effects 
along the line of sights should be estimated. Thirdly the noise due to the intrinsic asphericity 
of source galaxies should be treated properly. It is important to realize that the existence 
of noise can affect the detection of clusters considerably. Numerical studies (Hamana et al. 
2004; White et al. 2002) showed that the presence of noise reduces the efficiency of cluster 
detection significantly. Van Waerbeke (2000) investigated the properties of noise induced 
by the intrinsic ellipticities of source galaxies. He found that in the weak lensing regime, 
the lensing signal and the noise are largely uncorrelated if the smoothed convergence k is 
considered. Furthermore, to a very good approximation, the noise can be described as a 
two-dimensional Gaussian random field with the noise correlation introduced only through 
smoothing procedures. Then the technique of Bardeen et al. (Bardeen et al. 1986) can be 
used to calculate the number of peaks resulted purely from the noise. This provides us a 
possible way to estimate the contamination of noise on the abundance of lensing-detected 
clusters. The presence of noise also affects the height of peaks from real clusters. With 
numerical simulations, Hamana et al. (2004) tried to correct this effect empirically. In our 
future work, we will study the noise in great detail with the effort to establish a model to 
describe its effects on weak lensing cluster surveys. With the hope that this is achievable, 
we address in this paper the first two issues with the emphasis on the effects of the complex 
mass distribution of clusters themselves. 

Even for isolated clusters without any projection effect and without any noise, their lens- 
ing effects cannot be fully determined by their mass. Thus lensing-selected clusters cannot 
be truly mass-selected. Hamana et al. (2004) adopted the spherical Navarro-Frenk- White 
(NFW) (Navarro et al. 1996, 1997) density profile for a cluster to relate its smoothed peak 
K value with its total mass. Given a detection limit on k, they then obtained an ideal mass 
selection function with the redshift-dependent lower limit derived from the limit of k. The 
essence of their model is still that there is a one-to-one correspondence between the peak 
K of a cluster (at a given redshift) and its total mass. However, the complicated mass dis- 
tribution of a cluster can make its lensing effect much more complex than that predicted 
solely by its total mass. In this paper, we compare the peak k values of clusters from numer- 
ical simulations with the results by modeling the dark matter halos of clusters as isolated 
and spherical ones. Large differences in k exist for those clusters with the same expected 
K from the NFW-profile analysis. Our investigation further finds that the triaxiality of the 
mass distribution of clusters and large substructures in them contribute to this dispersion 
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significantly. 

The paper is organized as follows. In §2, we describe the lensing simulations and the 
cluster identification scheme from lensing convergence K-map. In §3, we present the results 
of our analyses. Summary and discussion are included in §4. 

2. Wecik lensing simulations 

In this paper, we consider the concordance cosmological model with = 0.3, Q\ = 0.7, 
r = 0.2 and as — 0.9, where Qrn ^-nd Q\ arc the present matter density and the energy 
density of the cosmological constant in units of the critical density of the universe, F is the 
shape parameter of the spectrum of the linear density fiuctuations, and ag is the rms of the 
extrapolated linear density perturbation smoothed over 8/i~^Mpc with h being the present 
Hubble constant in units of 100 kms~^ Mpc~^. 

The simulations we use (Jing & Suto 1998) have the box size of 100/i~^Mpc, and 256^ 
particles. The mass of each particle is 5 x lO^/i^^M©. The softening length is 39/i~^kpc. A 
typical cluster has a mass of ~ lO-'^^M© and a size of ~ Mpc, thus both the mass resolution 
and the force resolution are good enough for our studies. There are three simulation runs 
with different realizations of the initial conditions. For each simulation, it evolves from 
redshift z = 72 to z = with 1200 time steps equally spaced in terms of the cosmological 
scale factor. There are total 60 outputs with 14 of them between z = 1 and z = 0. Dark 
matter halos are identified with the FOF algorithm with the linking length b — 0.2 in units 
of the average separation of particles. 

In this paper, we focus on the weak lensing regime. In this limit, the shear 7 can be 
directly obtained from image distortions through e^^-* 7 + e^"^) where e is the ellipticity 
defined in the complex plane and / and S stand for image and source, respectively. Because 
both the convergence k and the shear 7 are determined by the lensing potential, k can be 
estimated from 7 (e.g.. Kaiser 1998). Then we have K,n{^) — i^i'^) + ^(^) where k is the true 
lensing convergence, and n is the noise part associated with the intrinsic ellipticity of source 
galaxies e*^*^^ (e.g.. Van Waerbeke 2000). Detailed analysis by Van Waerbeke (2000) shows 
that for the suitably smoothed k„, the lensing signal and the noise are uncorrelated in the 
weak lensing limit, and the noise can be well approximated by a two-dimensional Gaussian 
field. In his analysis. Van Waerbeke ignored the spatial clustering and the correlation of 
intrinsic ellipticities of source galaxies on arcmin scales. Therefore the smoothed noise is 
equivalent to the convolution of a point process by a smoothing window. According to the 
central limit theorem, the statistics of the smoothed noise is approaching Gaussian when the 
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number of source galaxies contained in the smoothing window is large enough. Although 

observations and numerical simulations indicate that the simplification mentioned above may 
be indeed valid on arcmin scales for source galaxies at z ~ 1 (e.g., Van Waerbeke 2000), the 
statistics of the noise can deviate from Gaussian if there are significant correlations in the 
spatial distribution or in the intrinsic ellipticity of source galaxies. In our future studies, we 
will investigate in detail the noise properties with these correlations included. 

It is very true that the existence of noise affects the identification of clusters significantly. 
On the other hand, for theoretical studies, it is fundamentally important to first understand 
the complexity of true lensing signals. The knowledge of this part forms the base for further 
investigations. Thus in this paper, we focus on true lensing signals with the emphasis on the 
influence of the complicated mass distribution of clusters and the projection effect. Therefore 
we mainly study noise-free k maps. The effects of noise on lensing cluster surveys will be 
carefully analyzed in our follow-up studies. 

In theoretical studies of lensing effects, generally one obtains the deflection angle a 
through ray-tracing simulations (e.g., Jain et al. 2000). Then the shear matrix is calculated 
by = dai/d(3j with (3 the source angular position. The convergence k and the shear 7 
can be derived from . In terms of denisty perturbations, to a very good approximation, 
we have (Jain et al. 2000) 



where g{x', x) = f^{x')f^{x~x') / f^ix) with r{x) being the comoving angular diameter distance, 
S is the density fluctuation at the actual position of photons, Hq is the Hubble constant and 
c is the speed of light. In the weak lensing limit, the Born approximation is valid and the 
integration in eq. (1) can be approximately performed along the unperturbed light path. 
We will see in §3.1 that our study presented in this paper concerns k values of order ~ 0.1 
for the Gaussian smoothing radius of 1 arcmin. Therefore the weak lensing approximation 
is well applicable. In this limit, k maps can be constructed from stacking different slices of 
mass distribution between sources and the observer together without intensive ray-tracing 
simulations. But we should keep in mind that at sub-smoothing scales, the k value can be 
much higher. The full lensing treatment should be carried out when studying lensing effects 
at these smaller scales. Our method of stacking in the weak lensing limit is detailed below. 

To generate a k map, we use a multiple-lens-plane algorithm. Many other similar algo- 
rithms choose line-of-sight directions to be parallel to the box edges (2;-direction, for example) 
(e.g.. White et al. 2002; Hamana et al. 2001). To diminish the effects of periodicity, the 
simulation boxes are randomly shifted when they are stacked together. In the algorithm we 
apply, the stacking is in a regular way without any shifts, but the line-of-sight directions are 
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chosen randomly. The effects of periodicity are minimal if we avoid special directions, such as 
30'', 45" and 60" (Jing 2004, private communications). Specifically, we fix the source redshift 
at Zs = 1. The observer is at z = and locates at the position {x, y, z) = (0.5, 0.5, 0.5). The 
comoving radial distance between z^ = 1 and z = is divided into 40 equal sections, whose 
corresponding redshifts z or the scale factors a can be obtained. Then each slice is filled by 
the outputting simulation boxes with their a value closest to the a value of the slice. The 
projected mass distribution in that shce is calculated on 1024 x 1024 grids for a 2 x 2 deg^ 
map. Thus along a line-of-sight direction, the convergence k is calculated through 



where i denotes different slices, 5i is the three-dimensional density perturbation, Oj is the 
scale factor, Aj is the comoving thickness of a slice, and are the comoving angular 
diameter distances to the source and to the lens, respectively, r^s is the comoving angular 
diameter distance between the lens and the source. 

We smooth k maps with a Gaussian window function of the form 



where 9g is the smoothing angular scale. Figure 1 shows a smoothed /t-map of 2 x 2 dcg^ 
with Oq = ^ arcmin. The grey-scale indicator shown in the right is in terms of the signal- 
to-noise ratio explained in the next section. The squares denote dark matter halos found 
in the corresponding directions with mass M > 10^^/i~^ Mq. Qualitatively, we do see good 
associations between peaks in the K-map and massive dark matter halos. 

To identify a match between a peak and a halo quantitatively, we adopt the algorithm 
of Hamana et al. (2004). Specifically, the matching is carried out in two directions: whether 
a peak has a corresponding halo and whether a halo has an associated peak in the K-map. 
Because it is not expected that the position of a peak coincides exactly with the central 
position of its corresponding halo, the matched pair candidate is searched within a radius 
of 2.88 arcmin around a peak or the center of a halo. The specific number 2.88 arcmin is 
chosen to be the same as that of Hamana et al. (2004), and it corresponds to 24 pixels in 
our K maps. The closest one is identified as the primary match. Then there arc five different 
peak-halo matching classes: double primary match, missing halo, false peak, secondary peak 
and secondary halo. A peak-halo pair is classified as a double primary match if both the 
peak and the halo are their corresponding primary matches. A missing halo is that there 
is no matched peak in the searching area around the halo. A false peak is that there is no 
associated halo around the peak. If a halo has a matched peak that has its own primary 
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matched halo, this halo is classified as a secondary halo. A secondary peak has an associated 
halo, but is not the primary peak of that halo (Hamana et al. 2004). Our investigation 
focuses on the class of double primary match. 

We do 2 X 2 deg^ surveys toward 720 directions for each simulation run. Of them, we 
discard those that have nearby halos {z < 0.1) more massive than 5 x 10^^h~^ Mq. The 
final number of surveys used in our analysis is 428, 487, and 401 for the three simulations, 
respectively. 



3. Results 



3.1. NFW density profiles 

The NFW density profile of a dark matter halo is described by 

Ps 



p{x) 



{r/rs){l + r/rs)2' 



(4) 



where and ps are the characteristic scale and density of the halo, respectively. Given the 
virial mass M^j^ of the halo, can be calculated through the fitting formula 



l + zllOi4/i-iM0/ 



(5) 



where the concentration parameter c„ = rvir/r-g with the virial radius of the halo, and 
c„* = 8 for the cosmological model considered here. Then ps can be obtained by 
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(6) 



Thus the density profile of a halo is fully determined by its mass. Cutting off the mass 

distribution of a halo at its virial radius, the surface mass density can be written as (e.g., 
Hamana et al. 2004) 

S(a;) = 2psrsf{x), 



(7) 



where x — r/vs and 
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(a;2-lKl+c„) ^ (l-a;2)3/ 
(x^-l)(l+c„) ^ (l-a;2)3/2d.lCCOb^(^_^^^^, 
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(x<l), 
(x=l), 

(1< X< Cn), 
(x> Cn). 
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The convergence n is then 

^^^^^sf{xl (9) 

where 

D 

^cr = — (10) 

with Dsi Dd and D^s are the angular diameter distances to the source, to the lens, and from 
the lens to the source, respectively, and 

Ks = 2p,rsKr- (11) 

Then one can obtain the smoothed central k from a dark matter halo by 

Ko^ J 2nK{9/9s)WG{9,9G)0d9. (12) 

Thus for a NFW spherical dark matter halo, there is a one-to-one correspondence between 
its Kq and its mass M^r- 

We consider double-primary-matched halo-peak pairs. For each halo, from its mass 
and redshift, we calculate the expected kqnfw, and compare it with the real k value of the 
corresponding peak. Similar to other studies (e.g.. White et al. 2002; Hamana et al. 2004), 
we use the signal-to-noise ratio u = Ko/anoise to represent the height of a peak with cr^oise 
given by 

9 ^ 1 

2 27r^|n,' ^ ' 

where cXe is the rms of the intrinsic ellipticity of source galaxies, Ug is the surface number 
density of source galaxies. Following Hamana et al. (2004), we take — 0.4 and Ug — 
30 arcmin"^. In most of our analyses, we use 9g = 1 arcmin. The corresponding cr„ojse = 
0.02. For a typical cluster of M ~ IO^'^^Mq, the peak k ~ 0.1 and thus z/ ~ 5 for = 
1 arcmin. Therefore it is appropriate to choose 9g in the arcmin scale for clusters of galaxies 
(Hamana et al. 2004). To illustrate the effect of using a different 9g, we will also present 
results with 9g — 2 arcmin, in which anoise — 0.005. 

In Figure 2, wc show the scatter plots of u^pw and Upeak for the three simulation runs, 
where i'nfw stands for the expected central signal-to-noise value for a halo modeled as a 
spherical NFW one, and Upeak is the halo's matched peak value from K-maps. We only 
consider halos with M > 5 x lO^^h"^ Mq and peaks with Up^ak > 3. From the plots, we do 
see a correlation between vnfw and Vpeaki but the dispersion is rather large. Quantitatively, 
we find that the dispersion a ~ 1.2, 1.3, and 1.4 for pnfw — 4.5, 5 and 6, respectively. We 
also note that the average vnfw ~ ^peak is slightly off zero, with i^mfw ~ ^peak ~ —0.3. 
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Hamana et al. (2004) found a positive bias with vnfw — ^peak ~ 0.24. It is likely that 
the small bias arises from simulation to simulation variances. On the other hand, the large 
dispersion indicates that for a cosmological model, the lensing-selected cluster abundance 
can deviate significantly from that predicted based on the spherical NFW description of 
cluster profiles. Therefore the NFW approximation can introduce large systematic errors in 
extracting cosmological parameters from lensing cluster surveys. 



3.2. TricLxial dark matter halos 

After a rough calculation, Hamana et al. (2004) attributed the dispersion to the sta- 
tistical distribution of the concentration parameter of the NFW density profile, and the 
projection effects of uncorrclatcd structures along the line of sight. Here we will study in 
detail the reasons of the dispersion, with particular attention paid to the triaxiality of dark 
matter halos. 

With high resolution simulations, Jing and Suto (2002) present a NFW-like triaxial 
density profile for dark matter halos, which is given by 

^^^^ ^ (i?/i?o)"a+^/^o)^-"' ^^^^ 

where pcr{,z) is the critical matter density of the universe at redshift z and 

«' = ^ + 7r + ^" (15) 

with the axial ratios Cx < Cy < 1. The characteristic scale i?o is related to the triaxial 
concentration parameter Cg = Re/ Rq with Re ~ OAbRmr defined in Jing and Suto (2002). 
Given the mass of a halo, the average Cx and Cy, and Cg, as well as their statistical distributions 
are derived explicitly from simulation results (Jing & Suto 2002). 

The effects of the triaxiality on strong lensing results, such as arc statistics and the 
probability of large-separation multiple images, have been studied in detail (e.g., Oguri et 
al. 2003; Oguri & Kccton 2004). The associated asphericity of the intraclustcr gas, revealed 
by X-ray and Sunyacv-Zcldovich effect observations, has also been investigated (e.g., Lee 
& Suto 2003, 2004; Wang & Fan 2004). Here we analyze the dispersion between i^peak and 
^NFW caused by the aspherical mass distribution of dark matter halos. 

The K profile of a triaxial dark matter halo can be analytically written as (Oguri et al. 
2003) 

<o = ^-^^^m. (16) 
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where the factor g is defined as 




with {6, 0) the polar coordinates of the hne-of-sight direction. The factor / is the same as 
that in equation (8) but with the variable ^ labeling the elliptical contours of the projected 
surface mass distribution. Specifically, 

e=-^{Ax' + Bxy + Cy'), (18) 

where we have assumed that the lensing plane is described by {x, y) coordinates, and 

A = cos^ e (\ sin^ 0+1 cos^ <p) + ^ sin^ ^> (19) 



5 = cos^sin20('^ - (20) 

^^x ^y' 



and 



C = ^ sin^ + 1- cos^ 0. (21) 

Thus for a triaxial halo of mass M^r at a fixed redshift its lensing signal depends on the 
line-of-sight direction. The statistical uncertainties of the axial ratios {cx and Cy) and of the 
concentration parameter (ce) also contribute to the dispersion of lensing signals. 

To demonstrate the dispersion caused by different viewing directions, in Figure 3 we 
show the simulation result (solid line) of the distribution of the associated peak v values of 
a halo of M = 3.9 x lQ^^h~^M.Q at z = 0.48 viewed along different line of sights. Also shown 
in the plot (dashed hne) is the smoothed central u distribution expected from a triaxial halo 
of this mass and redshift with its axial ratios determined from the real mass distribution 
of the halo. The concentration parameter Ce is taken to be the average value (Jing & Suto 
2002). It is seen that the lensing signal has a large dispersion, and the triaxial model does 
explain a large portion of it at high end. 

We now analyze the distribution of v expected from the triaxial model for a given lyNFW- 
Because vnfw depends both on the mass and on the redshift of a halo, halos of different 
masses at different redshifts can have same vnfw- Thus the conditional distribution is 

-..V / I V dz£^n[M^i,{iyNFW,z),z]p[iy\iyNFw-,M^ir{i^NFW:Z),z]dz 

Jo dz^n[M^;,[vNFW,z),z\ 
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where n(Mvir, z) is the mass function of dark matter halos. The probabihty function 
p[j^\i^NFw] M^irij^NFW, z), depends on the distributions of 9, 0, c^, Cy and Cg. Specifically, 

p[l'\l'NFW;Myi,.{l'NFW,z),z]dl' = ^jdCa:j d{Cx/Cy) j dCe j ^oi^^ ^ ^ 

]pice)pi0)pi<l>)}diy (23) 

where the distributions on c^, c^/cy, and Ce are taken from Jing and Suto (2002) [see also 
Oguri et al. (2003; 2004)]. The mass function of Jenkins et al. (2001) is adopted. 

In Figure 4, we show the u distribution for vnfw — 4.5, 5 and 6, respectively. The 
solid lines are the results from k, maps of three simulation runs. The dashed lines are the 
results of equation (22). The dash-dotted lines are the results for spherical halos taking 
into account the uncertainties of the concentration parameter. It is clearly seen that the 
dispersion solely from the distribution of the concentration parameter (dash-dotted lines) is 
small in comparison with the large dispersion shown in simulation data (solid lines). The 
triaxiality contributes additional dispersions, giving rise to a long tail at high i/ side. 

The high-end dispersion due to the triaxiality of dark matter halos has important effects 
on modeling the selection function for weak lensing cluster surveys. Since the lensing signal 
from a triaxial cluster depends not only on its mass and redshift but also on the viewing 
direction, the selection criteria corresponding to a detection limit Kum is no longer simply 
a mass limit Mii^{z) for each redshift z. Instead, for each Myj^, there is a probability 
p{Myir, z; Kiim) that its k. can be higher than Kum- Eq.(23) is the basic equation to obtain 
p(Mmr, z; Kiim)- Ouc of our futurc studies is to analyze p{Mmr, z; num) and further the weak 
lensing cluster abundances at different redshifts for different cosmologies. 

It is seen from Figure 4 that the triaxial model explains the dispersion better for more 
massive clusters. It is notably seen, however, that the range of u from simulations is broader 
than that predicted by the triaxial model, especially at the low end of u. For relatively weak 
detections {i^mfw = 4.5), the differences between the results from the model predictions 
and the simulations show up at both low and high u values. Possible reasons for these are 
discussed in the following. 

The differences between the simulation results and the triaxial model predictions can 
be attributed to two effects. One is the complex mass distribution of individual clusters that 
cannot be well approximated by the triaxial model. The other is the projection effect of 
uncorrelated structures along the hne of sight. To see the relative importance of the two, 
we generate lensing surveys for matched halos only. Specifically, only particles that belong 
to the matched halos are kept and all the other particles are removed from the simulation. 
In this way, we isolate the effect of the complexity of individual clusters. The lensing signal 
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from this artificial survey is referred to as Kgingie (or equivalently Vsmgie) ■ 

Figure 5 shows the scatter plots of fgingie vs. I'peak for ujqpw = 4.5, 5 and 6, respectively, 
where i^peak is the peak v value of matched halos from full simulations. Good correlations 
between the two are clearly seen. This association indicates that for these peaks, their k 
values are dominantly determined by the properties of their matched halos. The plots also 
reveal that the dispersion toward high i'peak is larger for smaller halos (lower i'nfw)- This 
shows, as expected, that the line-of-sight projection effects are relatively stronger for smaller 
halos. Quantitatively, we calculate the quantity < [{vpeak " ^single) / t^smgie]"^ >^^^ to estimate 
the strength of the projection effect. The values arc 0.25, 0.18 and 0.15 for u^fw = 4.5,5 
and 6, respectively. In the plots, the triangles denote clusters with more than one u peak 
associated with them. These clusters must contain large substructures. Note that we identify 
the u peak closest to the central position of a halo as its matched peak. For clusters with 
large substructures, the closest u peak is likely generated only by a substructure. Thus the 
vahic of matched i'peak underestimates the lensing signal from the full cluster. We indeed see 
that the triangles concentrate on the small i'peak part. This explains the lower-end extension 
of the distribution of i'peak in Figure 4. 

Now let us compare ^single with utH calculated from the triaxial model in which we 
determine the axial ratios and the orientation of a cluster from its real mass distribution. 
The results are shown in Figure 6. The symbols are the same as in Figure 5. First wc do 
see correlations between Vsingie and Vtrh with the associations better for single-peak clusters 
(crosses). It is noted, however, that the triaxial model predicts narrower scatters than Vsingie 
for relatively small clusters. For large clusters with pnfw = 6, the spread of vtH is about 
the same as that of ^single- Secondly, as expected, most of the triangles are in the lower right 
part of the plot. But there are several triangles scattered at the upper left part. We check 
individual cases, and find that multiple substructures along the line of sight can produce a 
higher lensing signal than that predicted by the smoothed triaxial mass distribution. 

We caution that the simulations we use are relatively small-sized (100/i~^ Mpc), and 
the number of large clusters is small. As we stack simulation boxes together to generate 
lensing maps, the same cluster can be surveyed many times along different directions. For 
i^NFW = 4.5,5, and 6, the number of matched peaks from lensing maps is, respectively, 
about 400, 300, and 200. However, the corresponding total number of clusters involved in 
our analyses is only 26, 23 and 14. The number of times that the most frequently surveyed 
clusters appear in our lensing maps can account for about 25% of the total number of matched 
peaks. Thus our statistical results can be influenced signiflcantly by the characteristics of 
only a few clusters. Similar analyses on larger simulations are desired for more robust 
statistical results. 
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We further do case-studies one by one. Of the respective 26, 23 and 14 clusters for 
J^NFW = 4.5, 5, and 6, the number of clusters with apparent substructures shown in lensing 
maps is 11, 12 and 10, with the corresponding fraction 42%, 52% and 71%. For these clusters, 
their lensing signals are mostly lower than those predicted by the triaxial model as explained 
earlier. A specific example is plotted in Figure 7. The 2;-projected mass distribution of this 
cluster is shown in Figure 8. Large subhalos are clearly seen. In fact, in this case, the cluster 
is in its early formation stage, and different substructures are just starting to merge. For 
those clusters without apparent substructures, most of them can be well described by the 
triaxial model. However, multiple hidden substructures along the line of sight can generate 
lensing signals higher than those from the smoothed triaxial model. An example is shown 
in Figure 9. For this cluster, the expected lensing signal from the spherical NFW mass 
distribution is i^new = 6. Taking into account the triaxiality of the mass distribution, 
the predicted lensing signals largely agree with that calculated directly from the real mass 
distribution of the cluster. For the extreme point where Vsingie is much higher than Vtri^ the 
mass distribution along this particular line of sight is shown in Figure 10. The upper panel 
is the real distribution and the lower one is from the smoothed triaxial model. Multiple 
major peaks are clearly seen in the upper panel, which explains the deviation of Vgingie from 
utri- Notice that when applying the triaxial model to individual clusters, we adopt the mean 
concentration parameter from Jing and Suto (2002) for a given set of axial ratios measured 
from the mass distribution of a cluster. The real concentration parameter can fluctuate 
around the mean. We do find that, in a few cases (5 or 6 out of 23 clusters with Vnfw = 5), 
the real mass distribution is much fiatter than that of the triaxial model with the mean 
concentration parameter, which results higher single than Vt^i- 

The above analyses were done for the smoothing radius Oq = ^ arcmin. To see the 
effects of we also perform studies with = 2 arcmin. Figure 11 shows the dispersions 
for 9g = 2 arcmin. Only the results from simulations (solid lines) and from the triaxial model 
(dashed lines) are shown. It is seen that the differences between the simulation results and 
the predictions of the triaxial model are larger than those for Oq — 1 arcmin. Quantitatively, 
for 9g — 1 arcmin and vnfw — 6, the dispersions are agimu ~ 1-5 and atn ~ 0.9, respectively. 
The corresponding numbers for Oq — 2 arcmin are agimu ~ 1-8 and atri ~ 0.6. In Figure 12, 
we show the scatter plots of Ugingie and Vpeak- First, as expected, the fraction of triangles is 
significantly less than that of = 1 arcmin. For 9g = 2 arcmin, the fraction is 9%, 19% 
and 11% for vj^pw = 4.5, 5, and 6, respectively. For 9g = ^ arcmin, the corresponding 
numbers are 23%, 32%, and 38%. The decrease of the effect of substructures does not 
however reduce significantly the dispersion of fpeak at low end for Oq = 2 arcmin. Secondly, 
the correlations of the two can still be seen in Figure 12, but the tightness of the association 
is worse than the case of = 1 arcmin. For Oq — 2 arcmin and vnew — 6, we have 
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Cr[{iysingle " J^peak) /j^single] ~ 0.24. For 9g ^ I arCHlin, a[{Usingle " J^peak) / J^single] ~ 0.15. The 

large dispersion reflects the severe projection contamination for Oq = '2 arcmin. In Figure 13, 
we show the particle distribution at several lens planes for a particular peak with Upeak = 4.5. 
This peak appears in both cases of ^g* = 1 arcmin and 6g = 2 arcmin. The corresponding 
primary cluster of the peak is located at plane 11. The three circles in each panel arc for 
angular scales of 1 arcmin, 2 arcmin and 3 arcmin, respectively. It is clearly seen that within 
the circle of 1 arcmin, the cluster itself is the dominant contributor to the lensing effect, and 
the projection effect from other planes is relatively weak. Within the circle of 2 arcmin, 
however, besides the cluster, the mass distribution on planes 7, 13, 31, 33 and 39 contributes 
significantly to the final lensing effect. Specifically, {vpeak ~ ^single) /^single ~ 0.37 and 0.8 for 
9g — i arcmin and 9g = 2 arcmin, respectively. We note that the projection effect causes 
large dispersions at both high and low ends. The joint effects of the stronger projection 
contamination and the reduction of substructures make the low-end dispersion not change 
greatly in comparison with that of = 1 arcmin. At the high end, however, the dispersion 
is significantly larger for 6g = 2 arcmin than for 6'g = 1 arcmin. Thus our analysis indicates 
that it seems too large to use 9g — 2 arcmin in weak lensing cluster surveys. 



4. Discussion 

Our study shows that for clusters of galaxies, their peak lensing signals appear much 
more complex than those predicted by the spherical mass distribution in which the total mass 
and the redshift of a cluster fully determine the strength of its lensing effects. Therefore the 
spherical model can introduce large errors to the prediction on the abundance of lensing- 
selected clusters. Those errors in turn will significantly bias the determination of cosmological 
parameters from future weak lensing cluster surveys. 

Hamana et al. (2004) attribute the complexity to the variation of the concentration 
parameter of the spherical NFW profile and the line-of-sight projection effect of uncorrelated 
structures. However, our analyses with 9g — i arcmin reveal that for relatively large clusters 
relevant to weak lensing cluster surveys, their lensing effects are mainly determined by the 
characteristics of individual clusters, and the line-of-sight projection effects play a minor 
role. For individual clusters, the triaxiality of their mass distribution is important and has 
to be taken into account in theoretical analyses on their expected lensing effects. In terms 
of the dispersion shown in Figure 4, the triaxiality contributes a considerable part of the 
high-end tail especially for massive clusters. The low-end extension is mainly attributed to 
the existence of large substructures in clusters of galaxies. Multiple substructures along a 
line of sight can also generate higher lensing signals than those predicted by the smoothed 
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mass distribution. Thus substructures also contribute, to a certain level, to the high-end 
dispersion. With = 2 arcmin, however, the projection effect is significant. Thus in terms 
of cluster detections with weak lensing surveys, the smoothing radius 9a — ^ arcmin is 
preferable to = 2 arcmin. 

Note that our analyses are based on the k maps generated under the Born approxima- 
tion. Therefore we may underestimate the projection effect because we integrate the weighted 
density fluctuations along the unperturbed light path. In our analysis with 6g = 1 arcmin, 
we find that the projection effect accounts for about 0.16 in the total dispersion of ~ 1.39 
of i^peak — J^NFW for i^MFW = 5. Hamana et al. (2004) quoted a value of ~ 0.5 from the 
projection effect estimated from the contribution of large-scale structures uncorrelated with 
halos (Hoekstra 2001). While the number 0.5 is not certain and the underestimate by the 
Born approximation in the weak lensing limit may not be as large as the numbers (0.16 vs. 
0.5) suggest, it is indeed desirable to quantitatively analyze the projection effects by full ray 
tracing simulations. 

In our analysis, we assume a fixed source redshift. For a real lensing survey, source galax- 
ies generally have a redshift distribution, which should be taken into account in modeling 
the lensing effects of clusters expected from the survey. 

With the mass distribution and the statistics for the triaxial dark matter halos from Jing 
and Suto (2002), the effect of the triaxiahty on lensing signals can be theoretically studied. 
Rather than simply mass selected, a much more complicated selection function for ^-limited 
lensing cluster surveys can be obtained. On the other hand, it is less straightforward to 
include substructures in theoretical analyses without numerical simulations. However, recent 
high resolution simulations have been able to present some statistics on substructures, e.g., 
their mass function and the spatial distribution (e.g. Gao ct al. 2004; Natarajan & Springel 
2004). These, in principle, allow us to analyze the infiuence of substructures on lensing 
effects without time-intensive simulations. 

Weak lensing effects are powerful and clean probes to the distribution of dark matter. 
Cosmological studies on lensing-selected clusters avoid physical processes related to the intr- 
acluster gas that complicate the cosmological applications of X-ray and SZ-selected clusters 
considerably. Our investigation reveals the complexities of lensing effects related to the mass 
distribution of clusters, which however, can be handled without fundamental difficulties. In 
a foreseeable future, weak lensing cluster surveys as well as cosmic shear observations will 
contribute greatly to cosmological studies. 

We sincerely thank the referee for the constructive and detailed comments and sugges- 
tions. We gratefully thank Y. P. Jing for kindly providing us simulation data and for useful 
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Fig. 1. — A 2" X 2° K-map from a simulation. 
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Fig. 2. — Scatter plots of Upeak vs. umfw foi^ the three simulation runs, respectively. 
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Fig. 3. — The distribution of u for a particular cluster with M 3.9 x lO^^/i~^M0 at redshift 
z — 0.48 {unfw — 6) viewed from different directions. The dashed hne is the Monte-Carlo 
result of the triaxial model for all the possible hne-of-sight directions. The axial ratios 
are calculated from the mass distribution of the cluster. The solid line is the result from 
simulations. 
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Fig. 4. — The distribution of v for all the clusters with Vf^p^ = 4.5, 5 and 6, respectively. 
The solid line is the result from n maps, the dashed line is the statistical result from the 
triaxial model, and the dash-dotted line is the result of shperical NFW model taking into 
account the statistical variation of the concentration parameter. 
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Fig. 5. — Scatter plots of Vpeak vs. Vsingie- The three plots are for vmfw — 4.5, 5 and 6, 
respectively. The triangles denote cases that multiple k single peaks appear within 2.88 arcmin 
around the center of the clusters. The crosses are for those that only one Kgingie peak exists 
within 2.88 arcmin around the center of the clusters. 




Fig. 6. — Scatter plots of vtH vs. fsmgie- The symbols are the same as in Fig5. 




Fig. 7. — The scatter plot of UtH vs. fgingie a particular cluster '95645' with unfw = 6- 
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Fig. 8. — The projected particle distribution along 2;-direction for cluster '95645'. 
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Fig. 9. — The scatter plot of Vtri vs. fgingie a particular cluster '71218' with unfw = 6- 
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Fig. 10. — The mass density along the hne-of-sight direction corresponding to the highest 
point in Fig.9. The upper panel is the real mass distribution of the cluster, and the lower 
panel is the one from the smoothed triaxial mass distribution. 
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Fig. 11. — Same as Figure 4 but with = 2 arcmin. Solid lines are results from simulations 
and dashed lines are the predictions from the triaxial model. 
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Fig. 12. — Same as Figure 5 but with 6g = arcmin. 
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Fig. 13. — The particle distribution in lens planes 7, 11, 13,31,33 and 39 for a particular 
peak with Upeak — 4.5. The three circles from inner most to the outer most correspond to 
angular scales of 1 arcmin, 2 arcmin and 3 arcmin, respectively. 



